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We show that the density matrix for a finite conductor attached to reservoirs obtained by Keldysh 
formalism is of MacLennan-Zubarev form. On the basis of the fact that the density matrix is the 
invariant part proposed by Zubarev, it is shown that Keldysh formalism may describe the irreversible 
processes and steady-state feature of the system. An important consequence of the MacLennan- 
Zubarev form of the density matrix is a generalization of the Kubo formula in a nonequilibrium 
case. On the basis of the result, we propose the formula of shot noise and a nonequilibrium identity 
between differential conductance, noise power and shot noise as a generalized Nyquist- Johnson 
relation. 
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1. INTRODUCTION 

Recently many-body problems in nonequilibrium 
steady states have been studied in mesoscopic devices. 
Keldysh formalismii£i£iii£ has been applied to study vari- 
ous systems. It has been used successfully to calculate ob- 
servable quantities such as differential conductance and 
noise power. 

However, in mesoscopic systems with electron corre- 
lation, progress has been limited even in the basic un- 
derstanding of nonequilibrium steady states. Of course 
the linear response theory is well established. In meso- 
scopic systems, the Kubo formula leads to the Nyquist- 
Johnson relation between conductance and noise power. 
However, there has been little discussion on the gener- 
alization of the Kubo formula and the relation between 
physical quantities in many-body systems in a nonequi- 
librium situation. 

One possible means of discussing general features un- 
der the nonequilibrium condition independent of differ- 
ent systems may be to consider density matrices. In this 
study, we consider a conductor in series with equilibrium 
reservoirs with a finite bias voltage. There is an approach 
to formulating the density matrix of this system&i, in 
which the density matrix is given by Keldysh formalism, 
although the term Keldysh formalism is not used. The 
obtained density matrix is similar to the equilibrium en- 
semble. The density matrix is expressed as an infinite 
series of operators, which are generated by the pertur- 
bative expansion of the s-matrix. Therefore, the density 
matrix has remained to be a formal expression. The gen- 
eralization of the Kubo formula under a nonequilibrium 
condition has not been successful. 

On the other hand, there has been progress by a dif- 
ferent method. By using a rigorous approach of the 
C* algebra, a density matrix of MacLennan-Zubarev 
form has been obtained for a finite system attached to 
reservoirs^. Concerning the transport properties of a 
one-dimensional noninteracting lattice or a quantum dot 
using the slave-boson technique, it has been reported that 



a current-current correlation function is given by the sum 
of differential conductance and an additional term. How- 
ever, the physical picture and origin of the additional 
term have not been clarified&i£. 

In this study we show that the Keldysh formalism 
of the density matrix of a finite conductor attached to 
reservoirs is of MacLennan-Zubarev form. We show that 
Keldysh formalism has sufficient features to describe the 
nonequilibrium steady state. We finally find that the 
MacLennan-Zubarev form of the density matrix leads to 
a nonequilibrium Kubo formula. On the basis of this 
form, we propose a nonequilibrium identity between dif- 
ferential conductance, noise power and shot noise, which 
is an extended Nyquist-Johnson relation in correlated 
systems in a nonequilibrium case. 



2. SYSTEM AND HAMILTONIAN 

We discuss a finite conductor attached to the left and 
right infinite reservoirs through boundary couplings. The 
system Hamiltonian is defined as 



H = H C + H L + H R + H, 



cL.R, 



(1) 



where H c , Ul,r and H c l.r describe the conductor, left 
and right reservoirs, and boundary couplings, respec- 
tively. The conductor part of the Hamiltonian is given 
by the sum of the noninteracting and interacting parts, 
H c = H c o + H c \. We stress that, in the following anal- 
ysis, the correlation effect is included. Generally, in the 
preparation of the Keldysh Green function, the system 
Hamiltonian H is divided into the nonperturbative and 
perturbative term a 11 ! 12 ' 13 as 



H = H + H 1 . 



(2) 



The nonperturbative term Hq = H c o+Hl+Hr describes 
three independent systems, which are the noninteracting 
conductor and infinite reservoirs. The perturbative term 
Hi = H c l^r + H c i is the sum of boundary couplings and 



the interaction inside the conductor. In many cases, H c \ 
is treated using a perturbation theory to calculate the 
Green function in Keldysh formalism. Notice that the 
following results do not depend on the method used to 
divide H into H q and Hi , which will be shown explicitly 
in section 4. 

Clearly, the energy scale of Hi is much smaller than 
that of Hq because Hi describes the effects in the finite 
conductor. Thus, it may be expected that the expansion 
of the s-matrix is a well-behaved series. 



3. ADIABATIC SWITCHING-ON OF H t IN 
KELDYSH FORMALISM 

Let us begin with the Keldysh formalism of a finite 
conductor connected to infinite reservoir s 11 ! 12 ' 13 . Gener- 
ally, the discussion based on Keldysh formalism is con- 
centrated on the Keldysh contour of the s-matrix and the 
Keldysh Green function. We carefully analyze the time 
dependence of Keldysh formalism in which the perturba- 
tive term is imposed adiabatically. 

In the initial state at to — — oo, the conductor and infi- 
nite reservoirs are decoupled in equilibrium state, where 
the initial density matrix is supposed to be given by 
po = e-^o-^c-^iVi-MRArR-no) and 0q ig the initial 

thermodynamic potential. The reservoirs are tuned to 
different chemical potentials. For simplicity, we define 

(pl + Ma)/ 2 = M = ° and (ml ~ Mfl)/ 2 = eV/2. We can 
rewrite the initial density matrix as 



Pa 



-l3(H -eV/2(N L -N R )-n a ) 



(3) 



In the following, we assume [Ho, Nl — Nr] — 0. 

The perturbative term is turned on adiabatically as 
ge~ c ' t 'Hi, where g is introduced to tune the magnitude 
of the coupling constant, and e is a positive infinitesimal 
number. We define the total Hamiltonian as 



H e = H + ge- e ^Hi. 



(4) 



The time evolution of the system obeys the Neumann 
equation 



i-p e (t) = [H e ,p e (t)l 



(5) 



where the initial state is given by eq. fl3J) , and [A, B] = 
AB — BA is the commutation relation. The formal solu- 
tion of the Neumann equation is given by 

p e (t) = U e (t,t )p e (to)U e (to,t) 

at S c {t,to)e lHoto p,{h)^ lHota Se{to,t)e lHot . 



,-iHni 



U £ (t,to) or its equivalent S e (t,to) satisfies 
i— U e (t,t ) = H e U e {t,t ), 



(6) 



(7) 



i-S £ (t,t ) = ge-^Hi(t)S e (t,t Q ), (8) 



where U e (t,t ) = e- iBot S e (t,t )e iHoto . In the fol- 
lowing analysis, we employ the explicit forms C/ £ (t,0) 
and S e (t, 0). Although the expressions of U e (t,to) and 
S e (t, to) are well known, here we show them explicitly for 
later use as 

US, to) = 1 + V i— f- / dti • • • dt n T[H e i ■ ■ ■ H en ], 

°° ( — \ n ) n f* 

5 e (t,i ) = l + V^ J ^ / dti---dt n 



e -<\^+-\'^T[Hi(ti)---Hi(t n ) 



(9) 



where T is the time-ordering operator and H en = Hq + 
ge-^\Hi. 

The expectation value of any operator O is given by 

Tr{Op e (t)} = Tr{G H (t)p e } = (0 H (t)), (10) 

where p £ = p e (0). The Heisenberg representation is de- 
fined by 

Q H {t) = U e (0,t)OU e (t,0) = S e (0,t)O(t)S e (t,0), 

(11) 

where 0(t) = g lH otQ e -iH t j g ^ e interaction represen- 
tation. The physical quantities are defined in the limit 
g — ► 1, e — > 0. On the basis of the observation that 
lim H e = H. In Keldysh formalism, O is supposed 

g— >l,e-*0 

to exist in the limit g — ► 1, e — > 0. 



4. MACLENNAN-ZUBAREV FORM 

Let us focus on the density matrix p t defined by p t = 
5 e (0, — oo)poS e (— oo,0), which may be rewritten as 

p t = cxp{-/35 £ (0, -Oo)(ff - eV/2(N L - N R ) 

-fi )5e(-OO,0)}. (12) 

This is simply a formal expression. In standard calcula- 
tions of Keldysh formalism, the s-matrix is almost always 
perturbatively expanded. In practice, we can sometimes 
estimate observable quantities based on some approxima- 
tions even if we do not have the complete knowledge of 
the density matrix. In this paper, we try to analyze the 
structure of the density matrix, which leads to a deeper 
understanding of the noncquilibrium properties of meso- 
scopic systems. 

Here, in order to study p e , we derive a key relation for 



A e = S e {0, -oo)AS e (-oo, 0), 



(13) 



where the operator A satisfies [A, Ho] = 0. We notice 
concerning A — Ho that a useful expression of A e is 



S e (0,-oo)H S e (-oo,0) 

. dS e (0,-oo) 



Ho + gHi - leg- 



9g 



S e (-oo,0), (14) 



which is obtained using the Gell-Mann and Low 
theore m 14 ' 15 . In the following, we derive a different 
expression for the arbitrary operator A that satisfies 



[A,H o ] = 0. 

First, we discuss the commutation relation between A 
and the s-matrix. 



J 



00 f-ini™ r° " 

[A, 5,(0, -00)] = Y, -^r~ / di i • • • d* n e- e (l tl l + - + l t -U(i 5 )- 1 nX;n-iCi-i 

n=l "■ " / -°° 1=1 

r[ffi(ti) • ■ • H l (ti- 1 )]J A (ti)T[H 1 (ti +1 ) ■ ■ ■ H x (t n )}, {h, • ■ • , ti-i} >U> {t l+1 , • ■ • , t n }. 



(15) 



r 



The "current" corresponding to the operator A is defined 
as 



J AH {i) 



d_ 
"di 



A H (t). 



(16) 



The calculation of eq. tfTS"]) can be performed as follows: 
We expand the s-matrix and consider the nth term in the 
series. Among the n\ possible time orderings, a certain 
term of t p > t q > • • • > t r is chosen. The commutation 
relation [A, H\{t p )H\{t q ) ■ ■ ■ Hi(t r )\ may be rewritten us- 
ing [i,ffi(t)j = (ig^JAit) for [A, Ho] = 0. With refer- 
ence to the time variable ti of J a, we can categorize times 
larger {t%, • • ■ , £z_i} and smaller {tj+i, • • ■ , t n } than t\. 
We sum possible time orderings except J ' Aiti) while tak- 
ing advantage of dummy integration variables. Finally, 
we obtain 



[A, 5 e (0,-OO)] = 



r° 

/ dte- e ^J A H{t)S e {0,~^)(17) 
J — 00 



It is easy to rewrite this as 



A £ = 5 e (0,-oo)A5 e (-oo,0) 

r° 

= A+ dte- eW J AH {t). 

J — 00 

Integration by parts in eq. (fl"8|) also gives, 

r° 

A e =e dte- eW A H (t). 

J — OO 



(18) 



(19) 



Note that A € is the invariant part of the operator A 
defined by Zubarev-i 6 -. The concept of the invariant part 
will play an essential role in the next section. 

By using eq. (fT5|) , the density matrix of eq. (fT2")) can 
be evaluated as 



expj-/3(# + / dte-^J eH (t) 



eV/2(N L -N R )-V I dte-^J cH {t) - fi ) J , 

(20) 



where J e H{t) and J c u{t) are the energy change and 
charge current, respectively, which are obtained by set- 
ting A = H and A = e/2(N L - N R ) in eq. (fTS]). 

The obtained density matrix is nothing but a type 
of MacLennan-Zubarev form of the density matri x 16 ' 17 , 
which is characterized by the time-integral of the "cur- 
rent" on the negative semi-infinite interval. We have ob- 
tained this form for the density matrix by dividing the 
system into Ho = H c q + Hl + Hr and Hi = H c l,r + H c i. 
On the other hand, for the case of Hq = Hl + Hr and 
Hi = H C Lji + H c i + H c q, it was obtained using a rigor- 
ous approach of a C* algebra^. Here, let us show that the 
obtained density matrix does not depend on the method 
used to divide the Hamiltonian. 

In eq. [fl8]). 



S,(0,-oo)H S e (~oo,0) = H Q + / dte' e ^J eH (t){21) 

J — 00 

seems to be dependent on the method used to divide H 
into Hq and H\. 

As pointed out previously, 5 C (0, — oo)HoS e (— 00, 0) is 
rewriten as eq. (fT4f . On the right-hand side of eq. (|14[) , 
Hq + gHi is clearly independent of the method used to 
divide H into Ho and Hi in the limit g — > l,e — > 0; 
hence, we discuss the last term 



Aa 



dS e (0,—oo) „ , „. ,„„. 

^9^H ^5 £ -oo,0. 22 

og 



Using 5 £ (0,-oo) = C/ £ (0, -00) and U e (0, -00) 
U e (— oo,0), we obtain 



ra=l 



(-i3) 



n— 1 poo 



/OO 
d*i---d* n T[# el •••#,„] (23) 
-00 



Now, we assume that H is given by the sum of the oper- 
ators A, B and C as H = A + B + C. We compare two 



cases: H 
From H en 



A4 
Ho 



B,Hi 



H fn = H„ 



■e\t„ 



C and H 



A, Hi = B + C. 



'Hi, we obtain 
(l-ge-^)B. 



(24) 



The last term (1 — ge e \ tn \)B which is dependent on the 
method used to divide the Hamiltonian vanishes in the 



limit g — ► 1, e — ► 0, hence J7 e remains the same. Thus, the 
density matrix in eq. (|20[) does not depend on a division 
of the Hamiltonian in the limit g — > 1, e — > 0. Therefore, 
we can conclude that the form obtained using the C* 
algebra is identical with the present form obtained using 
Keldysh formalism. 



5. CHARACTERISTIC FEATURE OF 
NONEQUILIBRIUM STEADY STATE 

The fact that p e is the invariant part of the density 
matrix leads to the understanding of the dissipation in 
the system and also of the mechanism needed to reach a 
steady state. These are the points clarified by Zubare\*i£. 
However, in view of its importance, we briefly repeat the 
essential parts of the arguments explicitly using Keldysh 
formalism. 



5.1. Invariant part of density matrix and causality 
condition 



The Neumann equation in eg. ([5]) is rewritten as 

i— lnp e (t) = [J3 r e,lnpe(<)]. 

By using eq.(|18j). lnp e can be also written as 



(25) 



lnp e = lnp + / dte-l £ l*J lnpojf (i). (26) 



where Jinp 0H (t) is defined using eq. (fT6|) for A = lnpo- 
Corresponding to cq. (fl"9|) . we rewrite it as 



die" |e|t In #,#(*) 



lim In p t = lim e 

e^O e->0 



= lim 

T a ->o 



1 



dthxp 0H (t). (27) 



We find that lnp e is given by the long-time average of 
the logarithm of evolution of the initial density matrix. 

This is nothing but Zubarev's definition of a nonequi- 
librium statistical operator—. He has pointed out that 
the entropy production is zero in the local equilibrium 
state; hence, the local equilibrium ensembles cannot de- 
scribe irreversible processes. To incorporate dissipation, 
the idea of coarse graining has been introduced. He has 
defined the noncquilibrium statistical operator as the 
long-time average of entropy which is logarithm of the 
local equilibrium ensemble. The time average has been 
defined at the negative semi-infinite interval, namely, the 
causality condition has been imposed. Thus, the entropy 
production has been proved to be positive, although the 
proof itself was limited in the linear response regime. 

Therefore, the adiabatic turning on of Hi precisely cor- 
responds to taking the invariant part with the causality 
condition. Therefore, we conclude that the adiabatic de- 
scription of Keldysh formalism may describe irreversible 
processes. 



5.2. Steady-state feature of time-correlation 
functions 



Based on the fact that p e is the invariant part, we will 
show that the steady state is described by the adiabatic 
switching on of Keldysh formalism in the limit e — > 0. 

We start our discussion with charge current. Follow- 
ing similar calculations to obtain eq. (|17[) . we show that 
the commutation relation between p e and t/ e (0,i) is ex- 
pressed as 



\p.Mo,t)] 



dtiJp sH (h)U e (0,t), (28) 



where Jp e n(ti) is defined by 
d 

Jp e H{tl) = —g-PeH(h] 



UcfrtJilp^HapeituO). 



(29) 



From eq. (|28p . we rewrite the expectation value of the 
current as 

Tv{p e J cH (t)} = Tr{U t (t,0)p e U € (0,t)J c } 

= Tr{p e J c }+ f At 1 TT{J- PcH {t 1 )J cH {t)}. (30) 
Jo 

We find that Jp £ i?(ii) is a key quantity for determining 
the time dependence of charge current. Jp l H{t\) is time 
derivative of the Heisenberg representation of the invari- 
ant part p e , which is one of the local integrals of motion 
defined by Zubarev-i£. 

Here, we discuss [p e ,H e i] in eq. ([2"9")) . By using the in- 
variant part of p e in eq. (|19|) . the commutation relation is 
rewritten as 

r° 

[p„H el ] = -ie / dt 2 e- e ^J PoH {t 2 ) 

J —OC 

+ e f dt 2 g e [p H(t 2 ),H €l ], (31) 

J — oo 

where we used 

[P0H(t2),H el ] 

= [poH(t 2 ),H e2 ] + [p 0H (t 2 ),H el - H e2 ] 

= -LWi 2 )+s(e- e|tl1 -e- e l t2 l)[poff(*2),ffi] (32) 

andg £ = 3e" c| * 2| (e" e|tl1 -e" e| ' 21 ). From eq . (fT8|) , the first 
term is given by 

- ie / di 2 e- e|t21 J PoH (t 2 ) = -Ufa - pa). (33) 

J — oo 

Due to g e , the domain of the integration of the second 
term in eq. (|31[) vanishes in the limit e — ► 0; hence, we can 
omit the second term. Thus, the commutation relation 
in the limit e — > becomes 



lim[p e , H el ] = lim -ie(p £ - p ). 

e— >0 e— >0 



(34) 



On the basis of the commutation relation, the second 
term of eg. ([50] ) is shown to be proportional to e, expressed 

as 

lhD.Tr{Jp m H(t 1 )J oH (t)} = limTr{i{p e , H el }J' cH (t u t)} 

£—►0 e— *0 

= lime-Tr{(p e - p )J' cH (*i,t)} 

- 0, (35) 

where J' cH (h,t) = U e (ti,0)J c H{t)U € (0,ti). The expecta- 
tion values of any operator for p e and po are assumed to 
exist in the limit e — > in Keldysh formalism. Therefore, 
e ■ Tr{- • ■ } vanishes in the limit e — ► 0. In other words, 
p e and H e i inside the trace commute with each other in 
the limit e — > 0. We have shown that hence, the second 
term of eq. (|3H| vanishes. 

Thus, we conclude that the expectation value of the 
current becomes independent of time in the limit e — > 0. 



\im(J cH (t)) = lim(J c ). 



(36) 



Following the same type of analysis, the time-dependent 
correlation function of charge current is obtained as 

lim(J cH (t)J cH (t')} = lim{J cH (t-t')J cH (0)).(37) 

e^0 e— »0 

The above two equations prove that the steady state 
is actually reached by the adiabatic switching-on of the 
perturbation term in Keldysh formalism. 



6. NONEQUILIBRIUM KUBO FORMULA 

The charge current and differential conductance are 
respectively defined by 

J(V) = lim (J cH (t)) = lim <J C >, (38) 

g— >l,e— >0 g— ► !,£— >0 



G(V) 



_ dJ(V) 

~ dV ' 



(39) 



In the last equality of eq. (|38|) , the steady-state feature of 
ea. ([5rJ|) is used. 

The differential conductance is obtained by differenti- 
ating the MacLennan-Zubarev form of the density matrix 
p e , eq. (|20p . with respect to V. 



G(V)= lim t3Ti\J c [e/2{N L -N R ) 

g-*l,e->0 



/V-.i, JsH(i) + ^) ft } 



(40) 

By substituting the identity 5 £ (0, — oo)S e (— oo, 0) = 1, 
dQo/dV is calculated as 

an 

-^L = -Tr{e/2(N L -N R )p } 

= -Tr{S e (0,-oo)e/2(N L - N R )S e (-oo,0)p } 

,0 

= -e/2(N L - N R ) - / dte- e|t| (Jcir(t)), (41) 



where in the last step eq. (|T8)) is used. Thus, we obtain 
the expression 



G(V) = lim /3 (5J c e/2(SN L - 5N R )) 

g-»l,e->0 



dte-^I^J^^W^W)), 



(42) 

where M = A - (A) and (ASB) = (JAJB) are used. 

Now, we rewrite the differential conductance into an- 
other expression. In preparation we derive several com- 
mutation relations. On the basis of [Hq, e/2(Nl — N r )] = 
0, 

S e (0, -oo)[Ho, e/2(N L ~ N R )}S t {~oo, 0) = (43) 

is obtained, which leads to 

H a + / dte~^dtJ eH {t), 

J — OO 

e/2(N L -N R )+ dte-^dtJ cH (t) 

J —OO 



(44) 



From this equation, one can derive 



f° 1 

p e ,H + / dte- eW dtJ eH (t) 

J — OO 

p e , e/2{N L -N R )+ dte-^dtJ cH (t) 

J — OO 



= 0. 



(45) 



Using the second commutation relation of eq. (|4"B")) for the 
correlation functions in eq. (I40|) , we can exchange current 
and the invariant part of e/2(A r L — N R ), which leads to 



G(V) = lim /3 [ (e/2(SN L - 5N R )SJ C ) 

g-*l,e-»0 



+ / dte^ tl (6J cH (t)SJ cH (0))) 

J — OO / 



(46) 

Moreover in cqs. (|42[) and (|4"6")l the time integrals can 
be rewritten at the positive semi-infinite interval with 
the steady-state property of eq. ([36|) . Thus, four distinct 
expressions of G(V) are obtained. Their sum, divided 
by four, gives the symmetrized expression of differential 
conductance, 

G(V)= lim P-({6J c ,e(SN L -6N R )}) 

g—>l,e—>0 4 

+ lim t ndte-^({5J cH (t),6J cH (0)}) 

S ->l,e-»0 4 J_ oc 

(47) 



where {A, B} = AB + BA represents the anticommuta- 
tion relation. 

We find that differential conductance is determined by 
two contributions: the current-current correlation func- 
tion and an unusual correlation function between current 
and the difference of the two particle numbers of reser- 
voirs. This result directly reflects the fact that the conju- 
gate of the bias V is the invariant part ofe/2(Nx,—N R ) in 
MacLennan-Zubarev form of the density matrix. Equa- 
tion (|47p is the main result of this study and may be 
considered as the generalization of the Kubo formula in 
a nonequilibrium situation for a general system of a finite 
conductor connected to the left and right reservoirs. 



7. PHYSICAL QUANTITIES 
7.1. Shot noise 

We introduce noise power as using 

/oo 
dte~^({6J cH (t),6J cH (0)})(48) 
-oo 

and define the current-charge correlation as 

S sh {V) = - lim ({SJ C , e(8N L - 5N R )}}. (49) 

3-*l,e->0 

By using S(V) and 5' s h(^), the differential conductance 
is given by 



g(v) = |s(v) - |s. h (n 



(50) 



Differential conductance and noise power are well-known 
observable quantities in mesoscopic systems. We propose 
that the 5 s h(V) given by eq. ([49|) be defined as shot noise 
in general, for both interacting and non-interacting sys- 
tems. 

To confirm that eq. (|l9")) is actually shot noise, we dis- 
cuss the simplest example of a noninteracting Anderson 
model for a dot attached to the left and right reservoirs, 



pka <r pka 



(51) 



where c v ^ a is the annihilation operator of the electron 
with the spin a in the left (p=L) or right (p=R) reservoir. 
d a is the annihilation operator at the dot and n a = d\d a . 
Vpka is the hybridization between the dot and leads with 
p = L,R. 

It is easy to calculate charge current using the Keldysh 
Green function 



/OO J 
—T(u)(f L -f R ), 



where f L<R = 1/(1 + e^ u ^ eV ^) are the Fermi distribu- 
tion functions. Transmission probability is defined as 



T(u) 



TrT 



Ll R 



2r 



r L + v R ( W - e d y + r2 ■ 



(53) 



where T LjR = ^2 k 2ir\VL,Rk<r\S(u - e L . Rk ) and T = 
(r^ + T R )/2. The differential conductance is obtained 
by differentiating J(V) with respect V as using 

/OO 1 
—T(u;)(f L (l~f L ) + f R (l-f R )). 

(54) 

On the basis of the definition of eq. (|4"5| . the noise 
power S(V) may be written as 

/>oo 1 

S(V) - 4e 2 / ^T(w)(/i(l - f L ) + /«(1 - /«)) 



2tt 



4e^ 



-T( W )(l-T( W ))(/i-/ R ) 2 . (55) 

Z7T 



In these way, we have confirmed that noise power is given 
by the sum of the equilibrium noise, the first line of 
eq. ([55|) . and shot noise. This result is well known in 
noninteracting system o 18 ' 19 . 

On the other hand, one can show that S^V) is given 
by 

/OO J 
5-T( W )(l-T(a;))(/ L -/ fl ) 2 (56) 
-oo 27T 

for a noninteracting system. This is nothing but the sec- 
ond term of eq. ([55|) . 

Therefore for the noninteracting system discussed here, 
we confirmed that 5 s h(V r ) expresses shot noise and that 
eq. ([50| actually holds. 



7.2. Generalization of Nyquist-Johnson relation 

In eg. ((50)) . we have found that conductance is given by 
the difference between noise power and shot noise. We 
would like to stress again that this general expression is 
also valid in interacting systems, and not limited to the 
noninteracting systems discussed in the literature so far. 
Therefore, we propose eq. ([50[) as a nonequilibrium iden- 
tity between physical quantities, which is a generalization 
of the Nyquist-Johnson relation into many-body systems 
in a nonequilibrium case. Below, we show that eq. (|50[) is 
actually a generalization of the Nyquist-Johnson relation. 

Let us first discuss the linear response regime. The 
density matrix at V — can be obtained using eq. (fi"4|) . 
It is well known that the last term of eq. (fj"4"|) describes the 
sum of all vacuum loops, effectively giving the correction 
of the thermodynamic potential Ail = Q — flo 1 ; hence, 
the equilibrium ensemble becomes 



(52) 



lim pJ v=0 = e-^ H -^. 

g->l,e->0 



(57) 



In any equilibrium ensemble, we can calculate the trace 
with the eigenfunctions of H\n) = E n \n) and J c = 
i[e(N L -N R ),H] so that 

S 8h (0) = -ie 2 ^(n|[(^iV L - 8N R ) 2 ,E n ]\n) = 0, (58) 

n=0 

is derived in general. S' s h(0) vanishes in the linear re- 
sponse regime. Thus, we obtain 



4fc B TG(0) = 5(0). 



(59) 



This corresponds to a Nyquist-Johnson relation in the 
linear response regime even for systems with electron 
correlatio n 18 ' 19 . 

This relation can also be derived using the fluctuation 
dissipation theorem. Thus the noncquilibrium identity 
cq. ([50|) shows that a naive generalization of Ak-QTG(V) — 
S(V) in the nonequilibrium case is not possible and there 
is the additional contribution 5 s h(V'), which plays an es- 
sential role in a nonequilibrium situation. 

Finally, let us discuss how to obtain shot noise experi- 
mentally. As we have discussed in the preceding subsec- 
tion for a noninteracting system, noise power may be di- 
vided into equilibrium noise and shot noise. Experimen- 
tally, we measure only the noise power. In actual cases, 
it has been discussed that noise power for eV > k^Tx 
includes shot noise significantly. On the other hand, the 
nonequilibrium identity suggests that shot noise as the 
difference between noise power and differential conduc- 
tance expressed as 



S ah (V) = S(V) - Ak B TG(V), 



(60) 



both of which are measurable quantities. Thus, we may 
use this identity as the definition of shot noise for both 
interacting and noninteracting systems. 



reservoirs is an example of MacLennan-Zubarev ensem- 
bles. On the basis of this fact, we pointed out that 
the adiabatic switching-on of the perturbative term of 
Keldysh formalism corresponds to taking the invariant 
part with the causality condition introduced by Zubarev. 
From this fact, we conclude that Keldysh formalism may 
describe an irreversible process and the steady-state fea- 
ture of the time-correlation functions For proof of this 
steady-state feature, it is essential to understand the fact 
that the density matrix and total Hamiltonian commute 
in the expectation values. 

Using the MacLennan-Zubarev form, we have derived 
the general expression of differential conductance, which 
is determined by the current-current correlation function 
and the correlation function between current and the dif- 
ference of the two particle numbers of reservoirs. We call 
it the nonequilibrium Kubo formula of a finite conductor 
attached to reservoirs. It is well known that noise power 
is given by the current-current correlation function. Wc 
have proposed that the correlation function between cur- 
rent and the difference of the two particle numbers of 
reservoirs gives the shot noise of the system, in general. 
This fact is confirmed by explicit calculations using the 
noninteracting Anderson model. Therefore, we propose 
that the nonequilibrium identity between differential con- 
ductance, noise power and shot noise is a generalization 
of the Nyquist-Johnson relation in general. 
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